Energy level statistics for models of coupled single-mode Bose-Einstein condensates 
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We study the distribution of energy level spacings in two models describing coupled single-mode 
Bose-Einstein condensates. Both models have a fixed number of degrees of freedom, which is small 
compared to the number of interaction parameters, and is independent of the dimensionality of 
the Hilbert space. We find that the distribution follows a universal Poisson form independent of 
the choice of couphng parameters, which is indicative of the integrability of both models. These 
results complement those for integrable lattice models where the number of degrees of freedom 
increases with increasing dimensionality of the Hilbert space. Finally, we also show that for one 
model the inclusion of an additional interaction which breaks the integrability leads to a non-Poisson 
distribution. 
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The application of random matrix theory to problems 
in physics can be traced back to the studies of Wigner 
[l|, y| in relation to the energy spectra of complex nu- 
clei. (For a recent survey of the current state of the sub- 
ject we refer to Q-) The mathematical formalism was 
largely provided by Dyson jj|, who determined that the 
energy level spacing distribution displays universal be- 
haviour depending only on the symmetry of the Hamil- 
tonian. (This is in contrast to the energy distribution, 
i.e. density of states, which is generically non- universal.) 

In this Letter we study the energy level spacing dis- 
tribution for two models of coupled single-mode Bose- 
Einstein condensates. Both of these models are integrable 
(i.e., there exists a set of mutually commuting operators, 
which includes the Hamiltonian, the number of which 
is equal to the number of degrees of freedom in the sys- 
tem). These models are also integrable in the sense of the 
Quantum Inverse Scattering Method (QISM) |^ based on 
the Yang-Baxter equation, as shown in [g. We find that 
the energy level spacing distribution follows the universal 
Poisson form P{s) = exp(— s), where s is the dimension- 
less energy gap parameter. Such a distribution indicates 
a lack of correlation between the energy levels, leading 
to random clustering. It was argued by Berry and Tabor 
[3, based on a semi-classical approach, that a Poisson 
distribution holds for integrable systems when the num- 
ber of degrees of freedom is greater than one. This result 
has been supported by many numerical studies on one- 
dimensional lattice models with many degrees of freedom, 
such as the Heisenberg model, the Hubbard model and 
the t — J model (at supersymmetric coupling) [a.l9l[l0l|. 
The integrability of these models results from the fact 



that each can be derived through the QISM. 

Our study for models of Bose-Einstein condensates is 
motivated by two factors. The first stems from the fact 
that the models we will analyse are integrable for a rel- 
atively large number of independent coupling parame- 
ters. This is in contrast to one-dimensional lattice models 
where integrability generally imposes severe constraints 
on the coupling parameters. The second motivation is 
that these examples are models with low numbers of de- 
grees of freedom acting in Hilbert spaces with arbitrarily 
large dimensions. In [l^ it is claimed that the reason for 
an integrable system to show a Poisson distribution for 
the level spacing is there exists a basis independent of the 
parameters in which the Hamiltonian is diagonal since 
there are as many commuting operators as the size of the 
Hilbert space. Or as explained in |lll | , if there are an in- 
finite number of conservation laws, then each subspace 
defined by a given set of quantum numbers contains a 
single level. Since the Poisson distribution occurs for the 
eigenvalue level spacings in random diagonal matrices, it 
should therefore apply to integrable systems. However, 
such an argument relies on there existing a large number 
of constants of the motion comparable to the dimension 
of the Hilbert space. The findings we present here sur- 
prisingly indicate that that the result is true even for 
systems with a small number of degrees of freedom, and 
does not rely on a semi-classical limit as in |3|. 

The first model we analyse is the two-site Bose- 
Hubbard model which has been widely applied for the 
study of two coupled Bose-Einstein condensates [IJ • The 
Hamiltonian is 

H = UnNi + UuNiN2 + U22Nl + ^iiNi+^i2N2 

C 

- -i^{a\a2 + alai), (1) 
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where the operators a^, a\^ Ni = a\ai, i = 1,2 are asso- 
ciated to two Heisenberg algebras with relations 



a^, a]] = 6ij, [oi, Oj] 
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TABLE I: Choices of coupling parameters for the Hamiltonian 
Q used to determine the density of states shown in Fig. 



The model describes Josephson tunneling between two 
condensates with tunneling strength £j/2, the parame- 
ters Uij are the amplitudes for S'-wave scattering and fii 
are chemical potentials. The Hilbert space of states is 
given by the infinite-dimensional Fock space spanned by 
the vectors 



Kn) = (alr(4n0) 



m, n = 0, 1,2, ...., oo. (2) 



As each basis vector is uniquely determined by the quan- 
tum numbers m, n, this model has two degrees of free- 
dom. The Hamiltonian commutes with the total particle 
number A^ = iVi -|- 7V2, so the existence of the two con- 
served quantities, H and N, shows that the model is 
integrable. 

In Ref. ^ it was demonstrated that (Q can be derived 
through the QISM and solved exactly using the algebraic 
Bethe ansatz. Here we will solve the model in the spirit of 
the co-ordinate Bethe ansatz, which allows us to exploit 
the tridiagonal structure of the Hamiltonian. Letting 
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it is easily shown, by directly evaluating the action of 
([3) on (|3J), that JVP) is an eigenstate of (Q provided the 
co-efficients am satisfy, for tti > 0, the recursion relation 
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such that ai = XgaQ, Qat+i 
£j{m+l)Xm _ 



and X„i is defined by 



iU_ 



12 



= UuN^ + (C/n + U22 - U12W 
2Uii)mN + fiiN + (/i2 - ^^l)m - E 



where E is the energy eigenvalue. It is clear that am 
is a polynomial in E of order m. The A^ + 1 roots of 
the equation a^+i = gives the complete energy spec- 
trum for the sector with N particles. Following 14] we 
can reduce the second order relation Q to a first order 
relation. Letting 
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with Yq 



1 and substituting into Q leads us to 
(m-iV-1) 1 
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FIG. 1: The density of states for the two-site Bose-Hubbard 
model with 1000 particles, and interaction couplings given by 
Table I. The approximately constant profile in the case of 
the couplings (a) lying in the Rabi regime is reflective of the 
semi-classical nature of the pendulum analogy given in |12 ]. 
In contrast, for the couplings (•) lying in the Fock regime the 
pendulum analogy is not semi-classical as indicated by a non- 
constant profile. The parameters (■) within the Josephson 
regime illustrate the crossover behaviour, where the density of 
states deviates from the semi-classical profile at low energies. 



showing that Ym admits a continued fraction expansion. 
This provides a convenient means to generate the poly- 
nomials am- For a given root E of aN+i = 0, this value 
may be substituted into the terms am, "m = 1,...,7V 
which permits us to compute the corresponding eigen- 
state through Q (although this is not necessary for our 
studies here). 

For U ^ Uii ^ U22 ^ —U12/2 it is common to divide 
the parameter space into three regimes [l2 |: viz. Rabi 
{U/£j « iV-i), Josephson (TV"! « U/£j « N) and 
Fock {N << U/£j). There is a correspondence between 
(^ in these limits and the motion of a pendulum |l2l |. 
For both the Fock and Josephson regimes the analogy 
corresponds to a pendulum with fixed length, while in 
the Rabi regime the length varies. For both the Rabi and 
Josephson regimes the system is semi-classical. By semi- 
classical we intend that the energy per particle forms a 
continuum in the limit of large particle number. The 
Fock case is not semi-classical (e.g., there is a finite gap 
in the energy per particle between the ground and first 
excited state; see llS],) and hence the argument of l2l| is 
not applicable. Fig. 1 shows the profile of the density of 
states n{E), supporting this picture, with the coupling 
parameters for the Rabi, Josephson and Fock regimes 
given in Table I. 

To study the full level statistics of this model in the 
most general context is prohibitive because of the infi- 
nite dimensionality of the Hilbert space. It is important 
to emphasize that one cannot simply restrict to a sec- 
tor of fixed particle number, as this has the effect that 
the restricted model has only one degree of freedom. In 
this instance, we have found that the Poisson distribu- 




FIG. 2: Energy level spacing distribution for the Bose- 
Hubbard model (0 in the sector 360 < N < 400, where 
the particle numbers are increased in units of 4. The 
total number of energy levels is 4,191. The distribution is 
independent of the choice of the coupling parameters, which 
are given by the first three columns of Table II. The inset 
shows the energy level spacing distribution for the remaining 
couplings given in Table II. The results show an excellent fit 
to the Poisson distribution P{s) = exp(— s), illustrated by 
the dashed curve in both the figure and inset. 
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TABLE II: Choices of coupling parameters for the Hamilto- 
nian Q used in Fig. |5| 



tion does not hold. As mentioned earlier, the Poisson 
distribution is only expected to hold when the numbers 
of degrees of freedom is greater than one. Therefore to 
investigate the level spacing statistics, we must conduct 
the analysis over a finite number of sectors with different 
particle number, in order to account for both degrees of 
freedom. 

We take subspaces of the Hilbert space comprised of 
sectors with fixed particle number N (each of dimension 
A^ + 1), starting at iV = 360, increasing in steps of 4 
particles, up to iV = 400. This gives a total number 
of 4,191 energy levels. We calculate these energy lev- 
els for the wide range of coupling parameters given in 
Table II. While the first three choices for the coupling 
parameters in Table II correspond to the Rabi, Joseph- 
son and Fock regimes, the remaining values were cho- 
sen randomly. Note that we have deliberately not taken 
Uii = U22 and fj,i = fi2 in all cases since this corresponds 
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TABLE III: Choices of coupling parameters for the Hamilto- 
nian © used in Fig. Inland Fig. 2] 



to a discrete symmetry upon interchange of labels 1 and 
2, leading to eigenvalue degeneracies which slightly com- 
plicates the analysis of the level spacing distribution. 

In determining the energy level spacing distribution 
for Fig. 2 no unfolding of the raw data was undertaken, 
nor required, unhke other studies [ESIlfl- ^^ ^^^h case 
the energy gaps were normalised by the largest gap, and 
a histogram built using, on average, 45 bins. In some 
cases there were a small number of large gaps (of ap- 
proximately two orders of magnitude larger) relative to 
the average gap, which were discarded. In all cases the 
number of discarded data points was less than 1% of the 
total. The curve y — 7exp(— /3s) was fitted to each of 
the data sets. Finally, each histogram was normalised 
by the factor 7"^ and the dimensionless gap parameter 
s was rescaled in each case by P~^. The data shown in 
Fig. 2 exhibits an excellent agreement with the theoret- 
ical Poisson distribution. 

The model we analyse next is one for an atomic- 
molecular Bose-Einstein condensate with two distinct 
species of atoms, denoted a and &, which can combine 
to produce a molecule c l.6| . It also has an interpretation 
as a model for second harmonic generation in quantum 
optics |15| where the non-linear terms in the number op- 
erators correspond to a Kerr non-linearity. The Hamil- 
tonian takes the form 

Ho = UaaN^ + UbbN^ + UccN^ 

+ UabNaNb + UacNaNc + UbcNbNc 

+tJ-aNa + fibNb + HcNc + n{a''b''c + c^ba) (6) 

which commutes with I = Na — Nb and the total atomic 
number J\f = Na + Nb + 2Nc. Along with Hq this es- 
tablishes that the model has three conserved integrals of 
motion. 

The model acts on the infinite-dimensional Fock space 
spanned by the vectors 



|;,m,n) = (at)'(6t)™(ct)"|0) 



(7) 



It is apparent that © has three degrees of freedom, spec- 
ified by the quantum numbers /, tti, n in |(7J), and is thus 
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FIG. 3: Energy level spacing distribution for the Hamiltonian 
@ in the sector with Af = 400 and J > 0, which has 20, 301 
energy levels. The coupling parameters are given by the first 
four columns of Table III, while the inset shows the data ob- 
tained for the remaining couplings in Table III. The distri- 
butions have universal behaviour independent of the choice 
of coupling parameters and follow the Poisson distribution 
P{s) = exp(— s), which is illustrated by the dashed curve in 
both the figure and the inset. 



FIG. 4: Energy level spacing distribution for the Hamiltonian 
||S|l in the sector 50 < TV < 100, where the particle numbers 
are increased in units of 10. This gives a total of 9,331 energy 
levels. The choice of coupling parameters is given by the 
first four columns of Table II, while the inset shows the level 
spacing distribution for the remaining coupling parameters. 
The results show a good fit to the Wigner surmise for a 
GOE distribution, which is indicated by the long-dashed curve 
(the short-dashed curve shows the Poisson distribution). 



integrable. For this model we can fix the total atomic 
number Af and make the restriction 2" > 0, which gives 
a Hamiltonian with two degrees of freedom acting in a 
finite Hilbert space of dimension d = {N"^ + 4Af + 3)/8 
for Af odd and d = {Af^ + 6Af + 8)/8 for Af even. 

The method described earlier for diagonahsing Q gen- 
eralises in a straightforward way to ©. Using this 
method we have diagonahsed © in the sector with 
Af = 400 and J > (this gives a total of 20,301 en- 
ergy levels) for the choice of couphngs in Table III, and 
determined the energy level spacing distribution in ex- 
actly the same manner as described for Q. The results 
depicted in Fig. O indicate a Poisson distribution for the 
level spacings which is independent of the choice of cou- 
pling parameters. 

Finally, we also examined the energy level spacings in 
the non-integrable Hamiltonian 



H = Ha + Hi 



(8) 



where Hi = a^fo+fc^a. Note that the presence of the term 
Hi means that Ti. commutes with Af, but does commute 
with 2. Each sector with fixed Af has dimension d — 
(A/'2 -f 4AA + 3)/4 for Af odd and d = {Af^ + 4Af + 4)/4 
for Af even. The results shown in Fig. 4 indicate that in 
this non-integrable case the level spacing distribution no 
longer fits the Poisson distribution, but lies closer to the 
Wigner surmise 
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which is the distribution for the Gaussian Orthogonal 
Ensemble (GOE) as expected for non-integrable time- 
reversal invariant Hamiltonians. A distinguishing feature 



of the GOE is level repulsion, which was well illustrated 
by Dyson 4\ who showed an analogy between the distri- 
butions for the level spacings of the GOE and a system 
of unit charges, with Coulomb repulsion, confined to the 
unit circle in two dimensions. 

In conclusion, we have shown that two integrable mod- 
els for systems of Bose-Einstein condensates exhibit uni- 
versal Poissonian behaviour for the distribution of energy 
level spacings independent of the coupling parameters of 
the Hamiltonians. Both models provide examples with 
a low number of degrees of freedom acting in an arbi- 
trarily large dimensional Hilbert space of states. To our 
knowledge this is the first stud y o f this type and com- 
plements previous studies d, y, [l3 of level spacing dis- 
tributions in integrable one-dimensional lattice models. 
For lattice models the number of degrees of freedom in- 
creases with increasing lattice length, and hence increas- 
ing dimensionality of the Hilbert space. Our results sup- 
port the view that all integrable quantum systems show 
a Poisson distribution for the energy level spacings, in- 
dicative of clustering due to randomness, provided the 
number of degrees of freedom is greater than one Q- 
Our analysis also shows that this result does not rely 
on the existence of a semi-classical limit, as it holds for 
cases where the non-linear interactions dominate, such 
as for the Fock regime of (QJ. We have also studied a 
non-integrable example and shown that for this case the 
Poisson distribution no longer holds, and level repulsion 
is displayed. However, it is important to stress that the 
converse is not true: examples of quantum systems not 
exhibiting level repulsion, which are nonetheless classi- 
cally non-integrable, have been shown to exist in [16|. 
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